Multi-phonon scattering and Ti-induced hydrogen dynamics in sodium alanate 
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We use ab initio methods and neutron inelastic scattering (NIS) to study the structure, ener- 
getics, and dynamics of pure and Ti-doped sodium alanate (NaAlH4), focusing on the possibility 
of substitutional Ti doping. The NIS spectrum is found to exhibit surprisingly strong and sharp 
two-phonon features. The calculations reveal that substitutional Ti doping is energetically possible. 
Ti prefers to substitute for Na and is a powerful hydrogen attractor that facilitates multiple Al-H 
bond breaking. Our results hint at new ways of improving the hydrogen dynamics and storage 
capacity of the alanates. 

PACS numbers: 61.12.-q, 63.20.Dj, 68.43.Bc, 81.05.Zx 



Developing safe, cost-effective, and practical means of 
storing hydrogen is crucial for the advancement of hy- 
drogen and fuel cell technologies. Presently, there are 
three generic routes for the solid-state storage of hydro- 
gen: (i) physisorption as in many porous carbon and ze- 
olite materials, (ii) chemisorption as in metal hydrides, 
and (iii) chemical reaction such as in complex metal hy- 
drides. Among the type (iii) materials, sodium alanate 
(NaAlH4) has received considerable attention because of 
its high hydrogen weight capacity and low cost. The 
release of hydrogen from NaAlH4 occurs via a two-step 
reaction: 



NaAlH4 
NasAlHe 



\ NasAlHe H 
3 NaH + Al 



I Al + H2 (3.7 wt%), 
- f H2 (1.8 wt%), (1) 



yielding a total of 5.5 wt% hydrogen. It was recently re- 
ported that a few percent of Ti doping in NaAlH4 renders 
accelerated and reversible hydrogen release under moder- 
ate conditions 1]. In spite of the extensive investigations 
of Ti-doped NaAlIl4 that have resulted, little is known 
about the mechanism by which Ti enhances the cycling 
kinetics of hydrogen [2, ly|- ^^di^ even the location 
of the Ti atoms remains unclear. While it is widely be- 
lieved that they reside on the surface of the material 1^, 
the possibility that Ti is substituted for Na has also been 
suggested js| , but convincing experimental or theoretical 
evidence is still lacking. 

Here we report neutron inelastic scattering (NIS) mea- 
surements of phonon density of states and first-principles 
total energy and dynamics calculations of pure and Ti- 
doped sodium alanates. We succeed in characterizing the 
main features in the observed spectra, which display sur- 
prisingly strong two-phonon contributions. Furthermore, 
the calculations show that it is most energetically favor- 
able for Ti to substitute for Na, breaking several Al-H 
bonds in its vicinity. We also find that Ti-doped NaAlH4 
can accommodate extra hydrogen near the dopant. In 
addition, we examine the effect of the Ti dopants on the 
vibrational spectrum of neighboring AIH4 groups, which 



could allow probing the Ti location in doped samples us- 
ing high resolution spectroscopic techniques. 

Three powder samples were investigated: pure 
NaAlH4, 2% Ti-doped NaAlH4, and NasAlHe- NaAlH4 
was prepared as described in Ref. 6. NaAlH4 was doped 
with 2 mol percent TiFa through mechanical milling ac- 
cording to standard procedure • Nas AlHe was synthe- 
sized and purified by the method of Huot et al. 8] The 
NIS measurements were performed using the Filter Ana- 
lyzer Neutron Spectrometer (FANS) at the NIST Center 
for Neutron Research 9] under conditions that provided 
energy resolutions of 2-4.5% in the range probed. 

The calculations were performed within the plane- wave 
implementation of the generalized gradient approxima- 
tion (GGA-PBE 0) to density functional theory (DFT) 
in the ABINIT package [llj. We used Troullier-Martins 
pseudopotentials treating the following electronic 

states as valence: 3(i and 4^ of Ti, 35 and 3p of Al, 3s 
of Na, and Is of H. We carefully tested the convergence 
of our calculations with respect to the plane-wave cut- 
off and k-point mesh. For example, for the Ti-doped 
alanate supercells (containing 96 atoms) we used a cut- 
off of 1250 eV and a 3x3x2 k-point mesh. The phonon 
spectrum of the pure compounds was computed using 
density functional perturbation theory (DFPT) 114] as 
implemented in ABINIT. We calculated the phonons cor- 
responding to two 2x2x2 q-point grids. An interpolation 
scheme was then used so that the powder- averaged NIS 
phonon spectra were computed for a 4 x 4 x 4 q-point 
grid within the incoherent approximation [l5| . 

For hydrogen dynamics studies, we first optimized the 
tetragonal NaAlH4 and monoclinic NasAlHe structures. 
The results, summarized in Tabled agree reasonably well 
with previous X-ray and neutron 0, studies. 
For example, the reported NaAlH4 lattice parameters 
(at 8 K) |li] are 4.98 and 11.15 A, respectively, while 
we obtain 4.98 and 11.05 A. The reported H position 
(0.237,0.384,0.547) is in good agreement with our result. 

We also calculated the energy change associated with 
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FIG. 1: Measured (top) and calculated (bottom) NIS spectra 
of NaAlH4. The calculated 1 and 1+2 phonon contributions 
are shown. The structure of NaAlH4 is shown in the inset; 
grey tetrahedra represent AIH4 units. 

the two reactions in Eq. (QJ and obtained, respectively, 24 
and 42 kJ per mol of H2 released, compared to the exper- 
imental enthalpies of reaction of 37 and 47 kJ/mol |l9j . 
The agreement is reasonably good taking into account (i) 
we cannot directly compare our energy differences with 
the experimental enthalpies (because we do not consider 
entropic contributions, etc.) and (ii) the error associ- 
ated with our first-principles approach is of the order of 
5-10 kJ/mol. This error is estimated from the error in 
the calculated cohesive energies of individual compounds; 
e.g.^ we get £^coh=4.53 eV for the H2 molecule, while the 
experimental value is 4.49 eV. 

The vibrational spectra of NaAlH4 and NasAlHe are 



TABLE I: Calculated structural parameters of NaAlH4 and 
NasAlHe. 
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FIG. 2: Measured (top) and calculated (bottom) NIS spectra 
of NasAlHe. The calculated 1 and 1+2 phonon contributions 
are shown. The structure of NasAlHe is shown in the inset; 
grey octahedra represent AlHe units. 

shown in Figs. Q and [3 respectively. The calculated one- 
phonon spectrum does not reproduce several experimen- 
tally observed features. Yet, including two-phonon pro- 
cesses yields results that are in excellent agreement with 
the observed spectrum of both compounds. Such strong 
multi-phonon contributions are unusual, but they seem 
to be typical of materials with AYRx groups. 

The nature of the different phonon bands in Fig. Q 
can be determined by computing the modes of the indi- 
vidual AIH4 groups embedded in the crystalline matrix. 
(Note that the NIS spectrum is dominated by hydrogen 
modes 15j.) As shown in Fig.IHb, we obtain four distinct 
groups of modes that correspond to the largest features 
in the one-phonon spectrum of Fig. Inspection of 
the eigenvectors allows the characterization of the modes. 
The lowest-frequency modes are AIH4 rotations, while 
the peak above 200 meV consists of stretching modes of 
the AIH4 tetrahedron. The modes in the two intermedi- 
ate sets are a mixture of displacements and stretches. 

Now we extend our calculations to investigate the pos- 
sibility of substitutional Ti doping of alanates. First we 
study whether Ti-doped alanate is energetically stable 
and, if so, where the Ti dopant goes. From the many 
substitutional and interstitial doping models that could 
be tried, we choose two that are experimentally moti- 
vated. Doping sodium alanate with solid TiCls by dry 
ball-milling results in formation of NaCl and partial des- 
orption of NaAlH4, which leads to formation of aluminum 
crystallites 20j. Hence, it seems pertinent to study the 
substitution of Al and Na by Ti. In the following we 
denote these doping models by "Ti^Al" and "Ti^Na" . 

We consider super cells containing 16 NaAlH4 formula 
units, and substitute only one of the Al or Na atoms by 
Ti. The cohesive energies are obtained as the sum of the 
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FIG. 3: Calculated one-phonon NIS spectra for (a) NaAlH4 
and for a single MH4 tetrahedron embedded in the lattice for 
various cases: (b) pure alanate, (c) Ti^Na, and (d) Ti^Al. 



individual atom energies minus the energy of the system. 
The results are given in TableHH E^oh of the pure alanate 
has been chosen as zero energy, so that positive entries 
indicate greater stability than the pure system. 

Note that the results in Table ^ give the energy 
change in reactions of the form Ti + Nai6Ali6H64 
Al + Nai6Ali5TiH64, which involve isolated atoms (Ti 
and Al in this case). These results thus measure the rela- 
tive stability of pure and doped systems. A positive (and 
large) entry in the Table indicates that, in principle, it is 
feasible to obtain the doped structure. Of course, in or- 
der to actually accomplish this, it is necessary to identify 
appropriate reactants that are able to release Ti easily 
and accommodate the atoms for which Ti substitutes. 

We find both Ti^Na and Ti^Al are energetically 
more stable than pure alanate; i.e. the system gains en- 
ergy by accepting a Ti dopant into the bulk and releasing 
a Na or Al atom. In addition, Ti^Na is found to be the 
most favorable substitution. 

The relaxed Ti^Na structure presents H atoms that 
come close to the Ti dopant. The shortest Ti-H distance 
is 2.05 A, to be compared with the 2.39 A Na-H distance 
in the pure system. Consequently, the distance between 
the hydrogens close to the dopant and their neighboring 
Al atoms is longer than the Al-H distance in the pure 
system; we obtain 1.70 and 1.64 A, respectively. This 
type of relaxation is to be expected, since Ti has more 
electrons than Na to share with neighboring hydrogens. 

It may seem surprising that Ti^Na has a higher co- 
hesive energy than Ti^Al. The typical valences of these 
atoms certainly suggest otherwise. However, Ti seems 
to be relatively large for the Al site. In the relaxed 
Ti^Al structure we get a Ti-H bond of 1.82 A, which is 
much longer than the 1.64 A Al-H bonds in pure sodium 
alanate. This size mismatch is the most likely cause for 



Ti^Al to be energetically less favorable. 

The above results do not change significantly when the 
electrons are allowed to spin polarize or when we allow 
both atoms and cell to relax (see E^^^^'^iSF) and E^"^^ 
columns in Table The spin-polarized calculations 

predict the Ti ion retains one unpaired electron in both 
Ti^Na and Ti^Al. 

A meaningful modification of the Ti^Na doping model 
is to introduce Na vacancies close to the Ti. Na vacan- 
cies should yield a more balanced sum of valence charges 
and have been invoked in the literature to argue that 
the Ti dopants reside in the bulk of the system [^. 
Our results for one and two Na vacancies are in Ta- 
ble IH denoted by "Ti^Na+Na^" and "Ti^Na+2Na^" , 
respectively. When spin polarization is allowed, the 
Ti^Na+Na^ structure is found to be more stable than 
the pure system, but significantly less stable than the 
doping models considered above. On the other hand, the 
Ti^Na+2Na^ structure is predicted to be quite unlikely. 

The dynamics of the neighboring H atoms could be 
used as a local probe for the Ti location. However, 
dynamical calculations for the whole 96-atom supercell 
are very computationally demanding. Instead, we cal- 
culated the vibrational spectrum of a single AIH4 group 
near the Ti dopant using the finite displacement tech- 
nique 21]. The results are shown in Fig. [2l In the 
Ti^Na case (panel c), the Ti dopant mainly affects the 
high-frequency modes, i.e. those involving stretching 
of the AIH4 tetrahedron. All the modes in that group 
soften. The mode that softens most, approximately from 
210 meV to 190 meV, is dominated by the displacement 
of H that is closest to the Ti, and essentially corresponds 
to its oscillation along the Al-Ti direction. This clearly 
indicates that the presence of Ti could facilitate breaking 
of the Al-H bond. 



TABLE II: Calculated cohesive energies, given in eV and per 
96-atom supercell. The result for pure NaAlH4 (231.922 eV) 
is taken as the zero of energy, ^coh™ is the cohesive en- 
ergy obtained by allowing the atoms to relax but imposing 
the NaAlH4 relaxed supercell; ^^^^^(SP) is the same but ob- 
tained from spin-polarized calculations; ^coh is the result ob- 
tained when allowing both atoms and cell to relax, and V the 
resulting volume of the 96-atom supercell in A^. 
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In the Ti^Al case (Fig. ISJi), the dynamics are mod- 
ified quite differently since we deal with a TiH4 group. 
This suggests that, by investigating the phonon spectrum 
of Ti-doped NaAlH4, one might determine whether Ti 
dopants go into the bulk of the system and, if so, where 
they are located. Motivated by this possibility, we mea- 
sured the phonon spectrum of a 2% Ti-doped sample, 
but obtained a result essentially identical to that of pure 
alanate shown in Fig. Q However, it should be noted 
that this does not rule out the possibility of substitu- 
tional doping in our sample since the amount of Ti is 
very small, and thus any dopant-induced feature in the 
spectrum should also be very small and hard to distin- 
guish from the noise. In addition, the NIS spectrum of 
pure sodium alanate presents significant two-phonon in- 
tensity in the 175-200 meV energy range (see Fig. 0, 
which makes it difficult to identify fine details. Higher 
resolution spectroscopic measurements, such as Raman 
scattering, might help elucidate this issue. 

The above results suggest that Ti dopants may facil- 
itate the breaking of the Al-H bond. We explored this 
possibility by moving one H atom to the immediate vicin- 
ity of the dopant and then relaxing the system. The re- 
sulting structure, which we denote by "Ti^Na (H)", is 
considerably more stable than the original Ti^Na dop- 
ing model (see Ecoh in Table Ej). In fact, we found that 
it has not one but two H atoms very close to the Ti. The 
shortest Ti-H distance is 1.81 A, and the corresponding 
Al-H distance is 1.89 A, i.e. 0.25 A longer than in pure 
NaAlH4. Thus the Ti dopant can indeed induce Al-H 
bond breaking, a necessary step for H2 release. 

These results indicate that the Ti^Na structure is a lo- 
cal minimum. In Fig. ^ we show the energy change along 
the transition path leading from Ti^Na to Ti^Na (H). 
The Ti^Na minimum is very shallow. The calculated 
energy barrier is around 2.5 meV~30 K, indicating that 
Al-H bonds would immediately break in the presence 
of Ti. Associated with the shallow well is a collective 
mode of the Ti^Na structure whose frequency can be 
roughly estimated to be ~80 meV. Yet, the frequency of 
the embedded AIH4 mode that we related to the Al-H 
bond breaking is about 190 meV (see the discussion of 
Fig-ISfc). The reduction from 190 to 80 meV is related to 
other atomic rearrangements, e.g., the second H coming 
close to the Ti atom, displacements of the AIH4 groups, 
etc. 

Along these lines, we also considered a less obvious 
possibility, namely, that Ti drags extra hydrogens into 
the system. Table El shows the result for "Ti-^Na+H", 
which corresponds to placing one extra H in the vicinity 
of the Ti atom. (In the calculation of this cohesive energy, 
the extra hydrogen was assigned one half of the energy of 
the H2 molecule, so that the resulting ^coh measures sta- 
bility against H2 release.) This structure turns out to be 
very stable. We find a Ti-H bonding distance of 1.82 A 
and several AIH4 groups approaching the Ti dopant. 
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FIG. 4: Energy along the path from Ti^Na to Ti^Na (H) 
structures (see text). Insets show the local structure and bond 
distances at the two minima. The zero of energy is arbitrary. 



In conclusion, we have used first-principles methods 
and neutron inelastic scattering to study pure and Ti- 
doped sodium alanate (NaAlH4), a material that holds 
great promise for reversible hydrogen storage. The total 
energy calculations indicate that substitutional Ti dop- 
ing in NaAlH4 is energetically stable. We find that the 
dopant prefers to substitute for Na and attracts several 
hydrogen atoms, softening and breaking the correspond- 
ing Al-H bonds. We also find it energetically favorable 
for the Ti to drag extra H atoms into the system. These 
results point to an interesting direction for future re- 
search, namely, the possibility of producing a new ma- 
terial, sodium-titanium alanate, that might benefit from 
the ability of Ti to accommodate extra hydrogens in its 
vicinity and thus exhibit improved H-storage capabilities. 

We acknowledge fruitful discussions with D. A. Neu- 
mann and Mei-Yin Chou. 
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